EXISTENCE OF DENSITIES FOR THE 3D NAVIER-STOKES 
EQUATIONS DRIVEN BY GAUSSIAN NOISE 

ARNAUD DEBUSSCHE AND MARCO ROMITO 

Abstract. We prove three results on the existence of densities for the laws 

of finite dimensional functionals of the solutions of the stochastic Navier-Stokes 
f^ I equations in dimension 3. In particular, under very mild assumptions on the 

O^l ■ noise, we prove that finite dimensional projections of the solutions have densi- 

^ I ties with respect to the Lebesgue measure which have some smoothness when 

measured in a Besov space. This is proved thanks to a new argument inspired 

by an idea introduced in [18] 
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1. Introduction 



r^ ! Our aim in this article is to give informations on the law of the solutions of the 

c^ \ stochastic Navier-Stokes equations in dimension three. The equations have the 

form 

on a bounded open set O. Here u is the velocity, p the pressure and v the viscosity 
Q \ of an incompressible fluid in the region O. The equations are supplemented with an 

-^ ' initial data and suitable boundary conditions. This equation has been the subject 

O ■ of intense researches, a survey can be found in the reference [11]. The forcing term 

^ \ fj in the problem above is a Gaussian noise which is white in time and depending 

on space (precise definitions will be given later). Under suitable assumptions on 
the noise, it is known that there exist weak solutions, both in the probabilistic and 
/\ , PDE sense. Their uniqueness is a completely open problem. Also, there exists a 

j^ \ unique strong solution on small time intervals. The situation is therefore similar 

to the theory of the deterministic Navier-Stokes equations. 

However, some more informations can be obtained if the noise is sufficiently non 
degenerate. It has been shown in [5, 9, 14, 15, 16] that it is possible to construct 
Markov solutions which depends continuously on the initial data. This indicates 
that the noise can be helpful to obtain more results in the stochastic case. 
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2 A. DEBUSSCHE AND M. ROMITO 

It is thus important to understand more deeply the imphcations of the addition 
of noise. In this article, we investigate the existence of densities of the law of the 
solutions. Existence of densities can be considered as a sort of smoothness of the 
solution, albeit a purely probabilistic one. 

A first classical difficulty is that the solutions live in an infinite dimensional 
space and that no standard reference measure, such as the Lebesgue measure in 
finite dimension, exists. It is tempting to use other measures as reference measure 
and, in [3, 6, 22] for instance, it is proved that for some equations the solutions 
have densities with respect to a Gaussian measure. Unfortunately, these results do 
not even cover the stochastic Navier-Stokes equations in dimension two. Another 
possibility is to try to prove existence of densities for finite dimensional functionals 
of the solutions. The problem of existence of densities of the solutions evaluated at 
a fixed spatial point has been already studied by several authors and we point to 
[23] for references. In the case of the two dimensional Navier-Stokes equation, finite 
dimensional projections of the solutions are studied in [21], where using Malliavin 
calculus, it is proved that there exist smooth densities. 

Unfortunately, it seems hopeless to use Malliavin calculus for the three dimen- 
sional Navier-Stokes equations. Indeed, it is not even possible to prove that the 
solutions are Malliavin differentiable. The reason for this is immediately apparent 
once one notices that the equation satisfied by the Malliavin derivative is essen- 
tially the linearisation of Navier-Stokes, and any piece of information on that 
equation could be used with much more proficiency for proving well-posedness. 

In this article, we propose three different approaches to the problem which give 
different results. First, in Section 3 we prove existence of densities under strong 
regularity and non-degeneracy assumptions on the noise (see Assumption 3.1). 
Due to the stronger assumptions, we are able to prove existence of densities for 
any smooth enough map of the solution with values in a finite dimensional space. 
In Section 4 we prove, by means of Girsanov's theorem, existence of densities for 
projections onto sub-spaces spanned by a finite number of Fourier modes, under 
the sole assumption that the covariance is injective (hence without regularity as- 
sumptions). A by-product of this technique is that the same statement holds also 
for the projection (onto the same sub-spaces) of the joint law of the solution eval- 
uated at a finite number of time instants (see Remark 4.4). Finally, in Section 5, 
under the same assumptions on the covariance of the Girsanov case, we prove 
again existence of densities with a completely different method, which extends an 
idea of [18]. This allows to show regularity of densities in the class of Besov spaces 
and to prove that they are in U spaces, for some p > 1, whereas the first two 
methods only provide existence of a density in L}. The regularity can be even 
slightly improved for stationary solutions. 

We believe that these result will helpful for future research on the Navier- 
Stokes equations. Recall for instance that the results of [21] have been a crucial 
step towards the fundamental result of [19]. Moreover, it seems that some ideas 
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introduced in this article are new and can be used to prove existence of densities 
in other situations, we refer to [8, 17] for two apphcations of the method used in 
Section 5. 

Acknowledgements. The second author wishes to thank the kind hospitality of 
the ENS de Cachan Bretagne, where part of this work was realised. 

2. Preliminaries 

We consider problem (1.1) with either periodic boundary conditions on the 
three-dimensional torus O = [0, 27i]'^ or Dirichlet boundary conditions on a smooth 
domain O C R'^. Let H be the closure in L^ = L'^{0; R^) of the space of smooth 
vector fields with divergence zero and satisfying the boundary conditions (either 
periodic or Dirichlet). The inner product in H is denoted by (■, ■) and its norm 
by II ■ \\h- The space V is the closure of the same space with respect to the 
H^ = H^{0; H?) norm. Denote by 11 the Leray projector, namely the orthogonal 
projector of L^ onto H. 

Let A = -nA, with domain D{A) = V f] if2(C,R^), the Stokes operator and 
let (Afc)fc>i and {ek)k>i be the eigenvalues and the corresponding orthonormal basis 
of eigenvectors of A. 

The bi-linear operator B : V x V ^ V is the (Leray) projection of the non- 
linearity {u ■'V)u onto divergence-free vector fields: 

B{u,v) = Il{u-Vv) , u,veV, 

and B{u) = B{u, u). The operator B can be easily extended to more general u, v. 
We recall that the following properties hold, 

(2.1) {ui,B{u2,U3)) = -{u-i,B{u2,Ui)) and {m, B{u2,Ui)) = 0, 

for all ui,U2, Ws such that the above expressions make sense. Moreover, there is a 
constant c > such that 

(2.2) \\A^B{ui,U2)\\h<c\\Aui\\h\\Au2\\h, Ui,U2eD{A). 

(see for instance [4]). We refer to Temam [31] for a detailed account of all the 
above definitions. 

We assume that the noise rj in (1.1) is of white noise type and can be described 
as follows. Consider a filtered probability space {fl, T , P, {J-t}t>o) and a cylindrical 
Wiener process W = Yliieisi Pi^ii where (/3j)igN is a family of independent Wiener 
processes adapted to (J-j)j>o and (gj)jgN is an orthonormal basis of H (see [7]). 
The noise f] is coloured in space with a covariance operator C G C{H) which is 
positive and symmetric. It is thus of the form fj = C^^. We assume that C is 
trace-class and we denote by a"^ = Tr(C) its trace. Finally, consider the sequence 
((T^)fcGN of eigenvalues of C. It is no loss of generality to assume that (gfc)fcgN is 
the orthonormal basis in H of eigenvectors of C: Cq^ = crlq;^. Further assumptions 
on C will be considered in the following sections. 
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With the above notations, we can recast problem (1.1) as an abstract stochastic 
equation 

(2.3) du + uAu + B{u) = C^ dW, 

supplemented with an initial condition u{0) = x E H . It is classical that for any 
X E H there exist a martingale solution of this equation. More precisely, there 
exists a filtered probability space (fi, J^, P, {Ji}j>o), a cylindrical Wiener process 
W and a process u with trajectories in C(R+; D(A-^))nL^,(R+, H)^rMf^^{K+ ; V) 
adapted to (J^)t>o such that the above equation is satisfied with W replacing W . 
See for instance [11] and the references therein for further details. 

The existence of martingale solutions is equivalent to the existence of a solution 
of the following martingale problem. We say that a probability measure P^ on 
C(R"''; D[A~^)) is a solution of the martingale problem associated to equation 
(2.3) with initial condition x E H ii 

.P,.[L-(R+,i7)nLL(R+;^)] = i, 

• for each (f> G D{A)), the process 

(6 - eo, 4>)+ I {^s, A^) - (B(6, 0), 6) ds 



is a continuous square integrable martingale with quadratic variation equal 

to t\\C'/'m, 
• the marginal of P^; at time is the Dirac mass 6^ at x, 

where in the formulae above {C,t)t>o is the canonical process on the path space 
C(R+;D(A"i)). 

The law of a martingale solutions is a solution of the martingale problem. Con- 
versely, given a solution of the martingale problem, it is not difficult to prove that 
the canonical process provides a martingale solution (see [11] for details). 

If K is an Hilbert space, we denote by C{K) the space of linear bounded op- 
erators from K into itself, hy Ttp : K -^ K the orthogonal projection of K onto 
a subspace F G K, and by span[xi, . . . ,x„] the subspace of K generated by its 
elements xi, . . . , Xn- Also B{K) is the set of Borel subsets of K. 

We shall use the symbol ^d to denote the Lebesgue measure on R'* and the 
symbol ^p to denote the Lebesgue measure on a finite dimensional space F in- 
duced by the representation by a basis. Finally, given a measure // and a mea- 
surable map /, we denote by /^yU the image measure of /i through /, namely 
{U^,){E) = ^{f-\E)). 

3. Existence of densities with non-degenerate noise: the 

Markovian case 

In this section we shall consider the following assumptions on the covariance. 

Assumption 3.1. There are e > and 6 G (1, |] such that 
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• Tt{A^+'C) < OO, 

• C-^A'^ e C{H). 

For example, C = A~" with a G (|, 3] satisfies the above assumptions. 

Under the above assumptions, it has been proved in [5, 9] and in [14, 15, 16] 
using a different method that there exists a family, indexed by the initial condition, 
of Markov solutions. 

We say that P(-, -, ■) : [0, oo) x D{A) x B{D{A)) -^ [0, 1] is a Markov kernel in 
D{A) of transition probabilities associated to equation (2.3) if P(-,-,r) is Borel 
measurable for every V e B{D{A)), P{t, x, ■) is a probability measure on B{D{A)) 
for every {t,x) G [0, oo) x D{A), the Chapman-Kolmogorov equation 

P{t + s,x,T)= I P{t,x,dy)Pis,y,T) 
Jd{a) 

holds for every t, s > 0, x G D{A), T G B{D{A)), and for every x G D(A) there is 
a solution F^: of the martingale problem associated to equation (2.3) with initial 
condition x such that P{t, x, T) = Pa;(6 ^ T) for all t > 0. 

Moreover, P{-,x, ■) and P^;, solution of the martingale problem, can be defined 
for all X G i^ and the Chapman-Kolmogorov equation holds almost everywhere in 
s. More precisely, for every x & H and every t > 0, there is a set / C R"*" such 
that the Chapman-Kolmogorov equation holds for all F G B{H). Also 

Vtip{x) = ^^-[v{^t)l t>0,xEH, 

defines a transition semigroup {Vt)t>o- It turns out that this transition semigroup 
has the strong Feller property, that is Vt(p is continuous on .0(^4) if (p is merely 
bounded measurable. Several other results can be found in the above references. 
In the following, by a Markov solution {Wx)x&h we mean a family of probabil- 
ity measures on C{IV~; D{A)) associated to transition probabilities satisfying the 
above properties. 

Let / : D{A) -^ R*^ be C^ and define, for our purposes, a singular point x of 
/ as a point where the range of Df{x) is a proper subspace of R*^. The following 
result is proved in Sections 3.1, 3.2 and 3.3 below. 

Theorem 3.2. Let Assumption 3.1 hold and let f : D[A) — )■ R*^ he a map such 
that the set of singular points (as defined above) is not dense. Given an arbitrary 
Markov solution {^x)x&h, l^t u{-,x) be a random field with distribution F^. Then 
for every t > and every x G H , the law of the random variable f{u{t] x)) has a 
density with respect to the Lebesgue measure =2^ on R^, which is almost everywhere 
positive. 

Example 3.3. We give a few significant examples for the previous theorem. 

• Functions such as f{u{t)) = \\u{t)\\'^jj, as well as any other norm which is 
well defined in D{A) admit a density with respect to the Lebesgue measure. 
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• In view of the results of the following sections, consider the case where / ~ 
TTp, where F is a finite dimensional subspace of D{A), iip is the projection 
onto F and / is given as f{x) = ((-, /i), . . . , {-Jd}), where /i, . . . , /d is a 
basis of F. Then the image measure (7ri?)#Pa; is absolutely continuous with 
respect to the Lebesgue measure of F. 

• Given points yi, . . . ,yd G R^ (or in the corresponding bounded domain in 
the Dirichlet boundary condition case), the map x H- {x{yi), . . . ,x{yii)), 
defined on D{A), clearly meets the assumptions of the previous theorem, 
and hence has a density on R'', since the elements of D{A) are continuous 
functions by Sobolev's embeddings. 

Remark 3.4. In view of [29, 1], the same result holds true under a slightly weaker 
assumption of non-degeneracy on the covariance of the driving noise. In few words, 
a finite number of components of the noise can be zero. 

By the results of [5, 9, 25] or [27], each Markov solution converges to its unique 
invariant measure. The following result is a straightforward consequence of the 
theorem above. 

Corollary 3.5. Under the same assumptions of Theorem 3.2, given a Markov 
solution (Px)xe-ff; denote by //^ its invariant measure. Then the image measure 
/#/x^ has a density with respect to the Lebesgue measure on R"^. 

Proof. If {P{t, X, ■))t>o,x£H is the corresponding Markov transition kernel and E C 
R*^ has Lebesgue measure J^d{E) = 0, then by Theorem 3.2 P(t, x, f~^E) = 
f#P{t, x,E) = for each x & H and t > 0. Then, by Chapman-Kolmogorov, 

f#^i.{E) = I P{t, X, f~^E) ^i^dx) = 0, 
Jd{a) 

since fi^{D{A)) = 1. D 

3.1. Reduction to the local smooth solution. Let x ^ C°°(R) be a function 
such that < X' < Ij xi^) = 1 for s < 1 and x(s) = for s > 2, and set for every 
R>0, XRis) = xii)- Set 

(3.1) Bniv) = XR{\\Av\\%)Biv) 
and denote by ur{-]x) the solution of 

(3.2) duR + [vAur + Br{ur)) dt = C^dW. 

with initial condition x G D{A). Existence, uniqueness as well as several regularity 
properties are proved in [16, Theorem 5.12]. We denote by Pr(-, ■, ■) and ¥^ the 
associated transition probabilities and laws of the solutions. 
Define 

(3.3) Tr = inf{t > : \\AuR{t)\\l > R}, 
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then again by [16, Theorem 5.12] it follows that t/j > with probability one if 
||Ax|||^ < R and that weak-strong uniqueness holds: every martingale solution of 
(2.3) starting at the same initial condition x coincides with Uf>{-]x) up to time t 
on the event {tr > t}, for every t > 0. 

Lemma 3.6. Let Assumption 3.1 he true and let f : D{A) — ?> R'^ be a measurable 
function. Assume that for every x G D{A), t > and R > 1 the image measure 
f^Pji[t, X, ■) of the transition density Pn^t, x, ■) corresponding to problem (3.2) 
is absolutely continuous with respect to the Lebesgue measure J^d on R"^. Then 
the probability measure f^P{t, x, ■) is absolutely continuous with respect to ^^ for 
every x E H every t > and every Markov solution {^x)x&h- 

Proof. Fix a Markov solution {^x)xeH and denote by P{t, x, ■) the associated tran- 
sition kernel. 

Step 1. We prove that each solution is concentrated on D{A) at every time 
t > 0, for every initial condition x in H. By [27, Lemma 3.7] 



E^ 



< oo, 



I ua'ds 

'-Jo -' 

for some 5 > 0. Thus P{s,x^D{A)) = 1 for almost every s G [0,t]. Recall that, 
for z G D{A) and r > 0, P{r, z, dIa)) = 1. We deduce that P{t - s, y, D{A)) = 
1, P{s,x,-)-a. s. for almost every s G [0,t]. Since the Chapman-Kolmogorov 
equation holds for almost every s, we have: 

P{t, X, D{A)) = IJ (^1 P{t- s, y, D{A)) P{s, x, dy)^ ds = 1. 

Step 2. Given x G D{A), s > and B C D{A) measurable, we prove the 
following formula, 

\Pis,x,B)-PR{s,x,B)\ < 2F4rR < s]. 

Indeed, by weak-strong uniqueness, 

P(s, X, B) = ¥F- [is eB,TR>s\+ ¥F- [is eB,Tn< s] 

= Pn{s,x,B)+E^^[is G B,rR < s]-E^^is G B,tr < s]. 

Hence the first side of the inequality holds. The other side follows in the same 
way. 

Step 3. We prove that the lemma holds if the initial condition is in D{A). Let 
B be such that ^d{B) = 0, hence PR{t,x, f-\B)) = for all t > 0, x G D{A) 
and -R > 1, then 

P{t + e, X, f'\B)) = j P(e, y, f~\B)) P{t, x, dy) 

<2 / P,[r«<e]P(t,x,d2/) + 2P(t,x,{||Ay||^>i?}), 

J{\\Ay\\„<R} 
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Since by [15, Proposition 11], F^Itj^ < s] i as s I ii \\Ax\\h < R, by first 
taking the limit as e J, and then as R^ oo, we deduce, using also the first step, 
that P{t + e,x, f~^{B)) — )■ as e J, 0. On the other hand, by [28, Lemma 3.1], 
P{t + e, X, f-\B)) ^ P{t, X, f-KB)), hence P(t, x, f-\B)) = 0. 

Step 4- We finally prove that the lemma holds with initial conditions in H. We 
know that P{t, x, f-\B)) = for alH > and x e D{A) if ^d{B) = 0. li x e H 
and s > is a time such that the a. s. Markov property holds, then 

P(t, X, r\B)) = j nt- s, y, r\B)) P(s, x, dy) = 0, 

since P(s, x, D{A)) = 1 by the first step. D 

3.2. Absolute continuity for the truncated problem. We now show that the 
law of f{uji{t] x)) has a density with respect to the Lebesgue measure on R'' for 
every R>0,xe D{A) and t > 0. We use Theorem 2.1.2 of [23] 

Let X G D{A), t > and R > 0. It is standard to prove that UR{t;x) has 
Malliavin derivatives and that V^UR{t] x) ■ qu = (Tkr]k{t, s; x), for s <t, where 77^ is 
the solution of 

(3.4) \ ( \ ^ ^ * 

Br is defined in (3.1), so that its derivative along a direction 9 is given as 

DBr{v)9 = Xr{\\Av\\%) {B{e, v) + B{v, 9)) + 2x'r{\\Av\\1){Av, A9)hB{v, v), 

and we recall that (g^, crk)ke'N is the system of eigenvectors and eigenvalues of the 
covariance C of the noise. Standard estimates imply that 

UR{t; x) G Di'2(Z}(A)) = |n : E[P^x||y + J^E / \\V'u ■ qkWl^^) < 00}. 

fc=0 "^° 

Moreover, for every /c G N and x G D{A), the function ^/^(t, s; x) is continuous in 
both variables s G [0, 00) and t G [s, 00). 

By the chain rule for Malliavin derivatives, the Malliavin matrix M.^ {t) of 
f{uR{t; x)) is then given by 

00 „( 
M{jit) = J2 {Dfi{uR{t] x))V'uR{t- x) ■ qk) {Df,{uR{t- x))V'uR{t- x) ■ qk) ds 

00 „i 

"X^^M {Dfi{uR{t;x))r]k{t;s,x)){Dfj{uR{t;x))r]k{t;s,x))ds 
k=i >^o 

for z, j = 1, . . . , 6?, where / = (J^, . . . , /^). 
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To show that A^-^(t) is invertible a. s., it is sufficient to show that if y G R*^ and 






ds 



k=l ■"" i=l 

is zero, then y = 0. This is clearly true, since if {M.^{t)y^y) = 0, then 

d 

y^ yiDfi{uR{t; x))r]k{t; s,x) = 0, P - a.s., 

i=\ 

for all /c G N and a. e. s < t. By continuity, the above equality holds for all s < t. 
In particular for s = t this yields 

d 

^yiDfi{uR{t;x))qk = 0, ¥ - a.s., 

i=l 

for all k > 1. Under our assumptions on the covariance, the support of the law of 
UFi{t]x) is the full space D{A) (this follows from Lemma C.2 and Lemma C.3 of 
[16]). Hence, uji{t]x) belongs to the set of non singular points of / with positive 
probability. We know that {qk)k>i is a basis of H, hence the family of vectors 
{Dfi{uR{t;x))qk,...,Dfd{uR(t;x))qk)k>i spans all R"* with positive probability, 
and in conclusion y = 0. 

3.3. Proof of Theorem 3.2. Fix an initial condition x & H and a time t > 0, 
and consider a finite-dimensional map / : D{A) — t- R'^ satisfying the assumptions 
of Theorem 3.2. Gathering the two previous sections, we know that f{u{t; x)) has 
a density with respect to the Lebesgue measure on R'^. 

Finally, to prove that the density is almost everywhere positive, it is sufficient 
to recall that under the assumptions of the theorem each Markov solution is irre- 
ducible (see for instance [10] or [16, Proposition 6.1]). 

4. Existence of densities with non-degenerate noise: Girsanov 

approach 

In this section we shall consider the following assumptions for the noise. 

Assumption 4.1. We assume that 

• the covariance C G C{H) is of trace-class, 

• ker(C) = {0}. 

We consider solutions of (1.1) obtained by Galerkin approximations. Given 
an integer A^ > 1, consider the sub-space Hn = span[ei, . . . , cn] and denote by 
T^N = t^Hn the projection onto Hn. It is standard (see for instance [11]) to verify 
that the problem 

(4.1) du^ + {uAu^ + B^{u^)) dt = ttnCUw, 
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where B {■) = 7iNB{TCj\f), admits a unique strong solution for every initial condi- 



tion x^ G Hm. Moreover, 



(4.2) Esuplln^ll^ 



<C,(1 + ||X^||^) 



for every p > 1 and T > 0, where Cp depends only on p, T and the trace a"^. 

li X & H, x^ = ttatx and P;^ is the distribution of the solution of the problem 
above with initial condition x^ , then any limit point of (P;^)Ar>i is a solution of 
the martingale problem associated to (1.1) with initial condition x. 

We prove the following result. 

Theorem 4.2. Fix an initial condition x & H and let F be a finite dimensional 
subspace of D{A) generated by the eigenvalues of A, namely F = span[e„j, . . . , Cnp] 
for some arbitrary indexes ni, . . . ,nF- Under Assumption 4-1 for every t > 
the projection Tipuit) has a density with respect to the Lebesgue measure on F , 
where u is any solution of (2.3) whose law is a limit point of the spectral Galerkin 
approximations defined above. Moreover the density is positive almost everywhere 
(with respect to the Lebesgue measure on F). 

Proof. Fix X E H and let w be a weak martingale solution of (2.3) with distribution 
Pa; and assume P;^*^ ^ P^;, where A/"/; '|" oo is a sequence of integers and for each k, 
P;^* is solution of (4.1). Fix a time t > and consider the Galerkin approximation 
(4.1) at level A^ > np. 

Step 1: the Girsanov density. We are going to use Girsanov's theorem in the 
version given in [20, Theorem 7.19]. Let w^ be the solution of 

dv"" + [uAv'' + S^(^^) - ^p5^(^^)) dt = TT^d dW, 

with the same initial condition as u^ . We notice in particular that the projection 
of v^ on F solves a linear equation (see (4.6) below) which is decoupled from 
f^ — TVpv^ . Moreover, it is easy to prove, with essentially the same methods that 
yield (4.2), that 



(4.3) E 



II A^IIP 

sup \\v 11^ 

iO,T] 



< CXD. 



Note that sup||^|| ^^^i{v,w) is a norm on tt^H, which is therefore equivalent to 
the norm of H on tt^H. We can then write: 

{■knB{v),w) = -{B{v,'Knw),v) < c\\w\\w^,o.\\v\\]j, 

therefore 

(4.4) \\tinB{v)\\h<cm\\v\\1, veH. 

Since the covariance C^ is invertible on ttnH, we deduce from (4.2), (4.3) that 

/ \\C~^TrNB{u^)\\jjds<oo, / \\C~^7rNB{v^)\\% ds < oo. ¥ - a.s. 

Jo Jo 
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By Theorem 7.19 of [20] the process 

(4.5) Gf = exp( / {C-hFB{u''),dWs) -\ f ||C-^7r^S(M^)||^ rfs 

is positive, finite P-a. s. and a martingale. Moreover, under the probabihty 
measure PAr(da;) = G^fN^du) the process 

Wt = Wt- I C~^ivfB{u^) ds 
Jo 

is a cyhndrical Wiener process on H and npu^it) has the same distribution as the 
solution z^ of the linear problem 

(4.6) dz^ + vAz^ = Tipd dW, 

which is independent of A^. In particular for every measurable E G F, 

FN[z^{t) eE]= FN[7rFu''{t) eE]= E^^ [Gf l£;(7r^M^(t))] 

and if ^f{E) = 0, where ^p is the Lebesgue measure on F, then lE{z^{t)) = 
0. Since Gf^ > 0, PAr-a. s., we have that lE{TTFU^{t)) = 0, PAr-a. s., that is 
FNliTpu'^^t) E E] = 0. In conclusion Tipu^it) has a density with respect to the 
Lebesgue measure on F. 

Step 2: passage to the limit. Now consider the weak martingale solution u of the 
infinite dimensional problem (2.3). We show that G^ '' is convergent. By possibly 
changing the underlying probability space and the driving Wiener process via the 
Skorokhod theorem, we can assume that there is a sequence of processes {u^'')k>i 
such that u^'' has law P^v^ and u^'' — )■ n a. s. in C{[0,T]] H^^) — where H^j is 
the space H with the weak topology — and in L^(0,T; H) for every T > 0. In 
particular the sequence {u^'')k>i is a- s. bounded in L°°{0,T] H) and thus a. s. 
strongly convergent in L^(0, T; H) for every T > and every p < oo. This ensures 
that G^ *" — 7- Gt, a. s., for every t, where Gt is the same as in (4.5) for u. 

Moreover, we notice that Gf > and finite a. s., since by (4.4) and (4.2), 

1 /"* 

- / \\C^ ^ Tip B {u)]]^^ ds < oo, a. s. 

2 Jo 
on the limit solution. 

Step 3: conclusion. We show that npu(t) has a density with respect to the 
Lebesgue measure on F. Let E G F with ^p{E) = 0, then for every open 
set J such that E G J we have by Fatou's lemma (notice that Ij is lower semi- 
continuous with respect to the weak convergence in H since J is finite dimensional) , 



E[GtlE{7Tpu{t))] < E[Gtlj{7rpu{t))] < 



< liminf E[Gf lj(7r^?/^(t))] = F[z^ (t) G J], 
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hence E[G'(lE(7rp'M(t))] = since J can have arbitrarily small measure and z^{t) 
has Gaussian density. Again we deduce that P[7rirM(t) E E] = from the fact that 
Gt > 0. Finally, the fact that the density of 7rpu{t) is positive follows from the 
results of [30]. D 

Remark 4.3. The bounds on the sequence (Gf^)7v>i are not strong enough to deduce 
a stronger convergence to Gt and hence to deduce the representation 

(4.7) E[cl){z^{t))] = E[Gt(P{7TFu{t))] 

in the limit, for smooth function : F — )■ R. Although this formula would provide 
a representation for the (unknown) density of T:pu{t) in terms of the (known) 
density of 2^, solution of (4.6), this would not characterise the law of irpu^t) 
by any means, since the factor Gt which appears in the formula depends on the 
sub-sequence {Nk)keN which ensures that Ptv^. ^ P- 
Vice versa, one could use the inverse density 

Gf = exp(- /" {C~"^7rFB{v''),dWs)-^f ||C-^7r^i?(t;^)||^ rfs 

which is also a martingale by [20, Theorem 7.19], to get in the limit 

(4.8) E[0(7r^n(t))]=E[Gi0(z^(t))] 

but the bound (4.3) for lA is not uniform in A^. 

Remark 4.4. Since our proof is based on Girsanov formula, it gives more infor- 
mation. In fact, it easily extends to show that for every ti,...,tm, the law of 
(iipu^ti), . . . ,7rpu{tm)) has a density with respect to the Lebesgue measure on 
F X ■■■ X F. 

Remark 4.5. In the proof of the above theorem we have actually used that the 
covariance is invertible only on F. Indeed Assumption (4.1) is too restrictive and 
the second property can be replaced by the following weaker property, 

• Fnker(C) = {0}, 

which allows a degenerate covariance. 

5. Existence of densities with non-degenerate noise: bounds in 

Besov spaces 

We now show that the density found in the previous theorem has a little bit 
more regularity than the one provided by the Radon-Nykodym theorem. At the 
same time we provide an alternative proof of existence of the density, which is 
based on an idea of [18]. We work again under the Assumption 4.1, although 
again Remark 4.5 may apply. 
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We prove in fact that the density belongs to a suitable Besov space. A general 
definition of Besov spaces i?p^(R'^) is given by means of Littlewood-Paley's de- 
composition. Here we use the equivalent definition given in [32, Theorem 2.5.12] 
or [33, Theorem 2.6.1] in terms of differences. Define 

{Alf){x) = f{x + h)-f{x), 
{Alf){x) = Al{Al-'f){x) = X;(-l)"-^(^)/(x + j/i) 

then the following norms, for s>0,l<p<oo,l<q'<oo. 



\\Ain\. dh^\ 



and for g = oo. 



A^ll 



B|,^ = ||;i|LP+SUp 

\h\<\ \l^ 



LP 



where n is any integer such that s < n, are equivalent norms of i?p^(R'^) for the 
given range of parameters. Note that if 1 < p < oo and s > is not an integer, 
then 5^p(R'^) = W'^^^K'^) (this is formula 2.2.2/(18) of [32]). We refer to [32, 33] 
for a general introduction on these spaces, for their properties and for further 
details on the topic. 

5.1. Besov regularity of the densities. We prove the following result. 

Theorem 5.1. Fix an initial condition x ^ H and let F he a finite dimensional 
subspace of D{A) generated by the eigenvectors of A, namely F = span[e„^, . . . , e„^] 
for some arbitrary indices rii, . . . , np. Under Assumption 4-i, for every t > the 
projection Tipuit) has an almost everywhere positive density fp^t with respect to the 
Lebesgue measure on F , where u is any solution of (2.3) which is limit point of 
the spectral Galerkin approximations (4.1). 

Moreover fp^t ^ Bl^^iTif-), hence fp^t G W^'^CR^), for every s G (0,1), and 
fp^t G L^CR'^) for any p G [1, ^), where d = dimF. 

Proof. Let u be a weak martingale solution of (2.3) with initial condition x and 
distribution F^, and assume P^* ^ Pa:, where A^^ t oo is a sequence of integers and 
for each k the probability measure P;^* is a weak martingale solution of (4.1) with 
initial condition ttn^x. Given a finite dimensional space F = span[e„j, . . . , e„^] 
and a time t, which without loss of generality is taken equal to t = 1, we wish to 
show that the random variable npu^l) has a density with respect to the Lebesgue 
measure on F (which we identify with R'', d = dimF). We first notice that, again 
by the results of [30], the density will be positive almost everywhere. 

For N > rip, let /at be the density of the random variable TVpu'^{l), where u^ 
is the solution of (4.1). The existence of /a? is easy to prove under Assumption 
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the indicator function of the interval 



4.1. For every e < 1, denote by r]^ = l[o,i 
[0, 1 — e]. Denote by m^'^ the solution of 

du^'' + {lyAu^'' + B{u^'') - (1 - r],)7rFB{u^'')) dt = n^d dW, 

where ttn is the projection onto span[ei, . . . , cn], and notice that u^'''{s) = u^{s) 
for s < t — e. Moreover for t G [1 — e,l], v = npu'^''^ satisfies 



(5.1) 



dv + uTTpAv dt = 11 pC 2 dW, 



Therefore, conditioned to J-*!, 
covariance 

Qf = 



e) = TTpu^'^l - e). 

7rp'M^'^(l) is a Gaussian random variable with 



Tipe 



VKpAs n unpA. 



Tip ds 



and mean 'Kpu^'''{l — e). We denote by g^^N its density with respect to the Lebesgue 
measure. Since Qp is bounded and invertible on F and its eigenvalues are all of 
order e, it is easy to see, by a simple change of variable (to get rid of the random 
mean and to extract the behaviour in e) and the smoothness of the Gaussian 
density, that 

\\ge,N\\Bl^ < Ce 2, 

holds almost surely with a deterministic constant c > 0. 
Fix (j) e C^iR'^), n > 1 and /i e R'^ with \h\ < 1, then 

E[{A-,<t>){npu^{l))] = n{K<P)i^Fu''{l)) - (A^Mnpu^'^il))] 

+ E[(A^0)(7r^n^'^(l))]. 

Consider the second term and use a discrete integration by parts 

E[(A^0)(^pn^'^(l))] = E[E[(A5:<^)(7r^«^'^(l))| J-i_, 



E 



E 



R'' 



Rd 



^h(P{x)ge,N{x) dx 

(t){x)/\\g^^N{x) dx 



<mLAhrn\9e,N\\Bi^ 

<c||0|U.oe-t||/i||". 



The first term of (5.2) can be estimated as follows 

|E[(A^c/>)(vr^n^(l)) - (A^0)(7rpn^'^(l))]| < 



ir^H ■:)E[(/)(7r^w 



j=0 

<c[0]«E[||7r^(«^(l) 



N 



.N,e, 



U 






<P{lTpu'''%l)+jh)] 
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where [0]q, is the Holder semi-norm of C^ilV^), and a G (0,1) will be suitably 
chosen later. Since 

7r^(n^(l) -n^'^(l)) = - [ e-"^(i-^)7ri.S(«^(s),«^(s))cis, 

Jl-e 

(4.4) and (4.2) yield 

E[|k^K(l)-n^'^(l))||] <cp r n\\u'^{s)\\j,]ds<CF{\\x\\l + l)e. 

Jl~e 

Gathering the estimates of the two terms gives 

|E[(A;:0)(7r^«^(l))]|<c[0],6° + ce-t||/i|n|0||^ = c||(/)||c;^| 



zn 

where the number c is independent of N, and we have chosen e = \\h\\ 2a+n. By a 
discrete integration by parts, 



n{A\c^){nFU^ {!))]= {A\c^){x)h{x)dx= {AlMix^x) dx 

jRd Jr.'' 

(we have switched from h to —h for simplicity) and so we have proved that for 
every h G R*^ with \h\ < 1, 

(5.3) / <P(y) ^^}[f^^ dy <c||0||c., 

with an = 2a "n " ^^ wish to deduce from the above inequality the following claim: 
The sequence {fN)N>np is bounded in Bj^ill^) for every 7 G (0, 1). 

Before proving the claim, we show how it immediately implies the statements of the 
theorem. Indeed, since BJ^{IV^) M- Lp(R'^) for every p G [1, -^] (by [32, formula 
2.2.2/(18)] and Sobolev's embeddings), it follows that the sequence {fNk)k>i is 
uniformly integrable and hence convergent to a positive function / G L^ilV^) 
which is the density of 7rpu{l) (we recall here that Ptvj. ^ P, where P is the law of 
u, hence the limit is unique along the subsequence {Nk)k>i). Moreover, since the 
bound in the claim is independent of A^, it follows that / G Bj^ilV^) for every 
7 G (0,1) and hence, using again the embedding of Besov spaces into Lebesgue 
spaces, / G Lp(R'^) for every p E [1, ^)- 

It remains to show the above claim. Let ip G iS(R'^), where iS(R'^) is the Schwartz 
space of smooth rapidly decreasing functions, and set cp = {I — Aaj'^^'^ip, where 
Ad is the Laplace operator on R*^ and [3 > a will be suitably chosen later. Notice 
that since ^"(R'^) = B^„^(R^) [32, Theorem 2.5.7, Remark 2.2.2/3] and since 
(/ — Arf)"'^/^ is a continuous operator from B^^{Il^) to B^ ooi^^) [32, Theorem 
2.3.8], it follows that 

C" < cU\\b^_^ < c||^|L<.-£ < Co 



■^00,00 
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where the last inequahty follows from the fact that L°°(R'^) M- B^^ilV^), since 
B^g(R'^) is the dual of 5f7"(R'^) [32, Theorem 2.11.2] and Sf7"(R'^) ^ L\'R'^) 
by definition, since (5 > a. 

Let gN = {I- ^d)-^'^fN. then (5.3) yields 

^lj{y){K9N){y) dy <c,\hr\mL-, 

hence ^Iqn G /.^(R'^) and 

Moreover, by [2, Theorem 10.1], Hfi-ArllLi < cWJnWl^ = c, hence {gN)N>np is a 
bounded sequence in S"^(R'^) and, since {I — A^)^^'^ maps Sf ^^(R'') continuously 

onto 5°^ (R'^) [32, Theorem 2.3.8], it follows that {fN)N>np is a bounded sequence 
in S";;^''(R'^) for every (3 > a. 

We notice that by suitably choosing n > 1, a G (0, 1) and f3 > a, the number 
a„ — /9 runs over all reals in (0,1): this can be easily seen by noticing that a„, — )■ 2a 
as n — )■ cxo. This proves the claim and consequently the whole theorem. D 

5.2. Additional regularity for stationary solutions. We can slightly improve 
the regularity of densities if we consider a special class of solutions, namely station- 
ary solutions. Consider again problem (4.1), it admits a unique invariant measure 
(see for instance [11], see also [26] for related results). Denote by Pat the law of 
the process started at the invariant measure. Every limit point is a stationary 
solution of (2.3), that is a probability measure which is invariant with respect to 
the forward time-shift (other methods can be used to show existence of stationary 
solutions, see for instance [12]). 

The idea that stationary solutions may have better regularity properties has 
been already exploited [13, 24]. 

Theorem 5.2. Let F be a finite dimensional subspace of D{A) generated by 
the eigenvalues of A, namely F = span[en^, . . . , 6^^] for some arbitrary indices 
ni, . . . ,np. Let u be a stationary solution of (2.3) which is a limit point of a se- 
quence of stationary solutions of the spectral Galerkin approximation. Under As- 
sumption 4-1, the projection npull) has a density fp with respect to the Lebesgue 
measure on F , which is almost everywhere positive. 

Moreover fp € Bl^i^'^), which in particular implies that fp G W^'^i^'^) for 
every s G (0, 2), where d = dimF. 

Proof. We proceed as in the proof of Theorem 5.1. Fix a stationary solution u with 
law P, a sequence Pat^. ^ P of stationary solutions of (4.1) and a finite dimensional 
space F = span[e„j, . . . , e^^]. Write again 

E[(A^'0)(vr^«^(l))] = E[iAl<j>)irrpu^il)) - (Ai:c/>)(vr^n^'^(l))] 

+ E[(A5:0)(vr^n^'^(l))]. 
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where this time w^''^ is defined as the solution of 

+ (1 - 7],)TTFB{e^^^'^^+'^ u^''(l - e))) ds = tin^ dW^ 

so that again u^{t) = u^''^{t) for t < 1 — e, and for t > 1 — e the process Tipu^''^ 
satisfies 

dv + [v-kfAv + 7iFB{e-'^^'-^+'^ u^''{l - e))) ds = upd dWs, 

which is the same equation as in (5.1) with an additional adapted external forcing. 
As before, conditioned to -Fi-e, 'Kfu^''^{^) is Gaussian with covariance Qp- Thus, 
the second term of (5.4) has the estimate 

|E[(A^0)(7r^^^'^(l))]|<c6-t|/,|«||0|U. 

We claim that 

(5.5) ^[hFu'^iX) - ^f«^'^(1)||h] < cei. 

Before proving (5.5), we show how to use it to conclude the proof. Indeed, as 
before, the first term on the right-hand side of (5.4) is bounded from above as 

|E[(A^0)(W^(1)) - (A^0)(w^'^(l))]| < c[0]„E[|| w^(l) - T^FU^'^^mHT 

and so 

f (A;:/;v)(x)0(x)rfa:<c[0],et" + ce-t|/in|0|U<c||</)||c.|/i|"^ 

by choosing e = |/i|2"/(3a+n)^ where this time a„ = ^^. As in the proof of 
Theorem 5.1, the above estimate yields that {fN)N>nF is a bounded sequence in 
Bf ^, for every s < a„ — a. Since a„ — a — )■ 2q: as n — >■ oo and a G (0, 1) can be 
arbitrarily chosen, we conclude that {fN)N>nF is bounded in Bf^ for every s < 2. 
In particular, since Bf^^CR^) M- Bf-^^CR^) = W^'^^R^), {fN)N>nF is also bounded 
in W^'^CR^) for every s < 2. 

We conclude with the proof of (5.5). We have that 



1 

l-e 



7r^(ix^(l)-n^'^) = / e-'^^'"'KF{B{e-'^^'-'+'K''{l-e))-B{u''{s)))ds, 
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hence by (4.4) and Holder's inequality, 
(5.6) 

E[|ki.(«^(l)-n^'11))||] 



< 



< c / E[(|| e-^^(^-i+^) «^(1 - 6)11^ + ||n^(.)||^)|| e-^^(-i+^) n^(l - e) - n^(.)||^] ds 
n-e 



<c I E[{\\u''{l-e)\\H + \\u''{s)\\H)TE 



.^A(s-l+e) ^N 



4 , 3 



M^(l-e)-«^(s)|||]^ds 



l-e 
1 



<c / E\\\e-'^^'-'+'^u''(l-e)- 



u'^is) 



l-e 



\lVds, 



since E[||w^(s)|||^] is finite, constant in s and uniformly bounded in A^. Now, for 
5G(l-e,l), 



-uAis-l+e) ^N 



N I 



?/" l-e -m" s 



Jl-e Jl-E 

and so 

E[|| e-'^^(^"i+^) w^(l - e) - n^(s)|||] < 



< E 



|e-"^(^-'')S(u^(r))|Lrfr 



l-e 



E 



^_,A(,_,)^i^ 



l-e 



f/. 



= S + @- 

To estimate [i] we use the inequality 

\\A-'^Biv)\\H < c\\v\\% < c\\v\\i\\v\\l., 

standard estimates on analytic semigroups and we exploit the fact that u^ is 
stationary. 



0<E 



< cem 



'^ " .N/„\\\h ]\„.N, 



1 



\u^ir)\\l\\u'\r)\\ldr 



i-€ [s — r)3 

.N\\l |L,Af||2 



u^{r)\\Mu''^ir)rvdr 



cesElWw WhWu 



2 

ce3. 



h\\un\v] 



The second term is standard, 



§<E 
and in conclusion 



^_,^(,_,)^i 



l-e 



H. 



<(^eTr(C))^=c6§, 



E lie 



-vA{s-l+e) ^N 



u^{l-e)-u^{s)\\l]<ce-^, 
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hence from (5.6), 

which proves (5.5). D 
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